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Abstract: For over a decade, the United States Department of Energy, and engineers, 
geologists, and scientists from all over the world have investigated the potential for 
reducing atmospheric carbon emissions through carbon sequestration. Numerous reports 
exist analyzing the potential for sequestering carbon dioxide at various sites around the 
globe, but none have identified the potential for a statewide system in Florida, USA. In 
2005, 83% of Florida's electrical energy was produced by natural gas, coal, or oil 
(e.g., fossil fuels), from power plants spread across the state. In addition, only limited 
research has been completed on evaluating optimal pipeline transportation networks to 
centralized carbon dioxide repositories. This paper describes the feasibility and preliminary 
locations for an optimal centralized Florida-wide carbon sequestration repository. Linear 
programming optimization modeling is used to plan and route an idealized pipeline 
network to existing Florida power plants. Further analysis of the subsurface geology in 
these general locations will provide insight into the suitability of the subsurface conditions 
and the available capacity for carbon sequestration at selected possible repository sites. The 
identification of the most favorable site(s) is also presented. 

Keywords: CCS; carbon sequestration; centralized repository; transportation optimization; 
Florida; storage alternatives 
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1. Introduction 

The potential for global climate change due to warming of the planet has been studied and debated 
over the last decade. The majority opinion is that the climate is changing due to emissions of carbon 
dioxide and other greenhouse gases. Scientists from around the globe generally agree that reducing 
greenhouse gas emissions is of paramount importance. One way to achieve sizable reductions over the 
next few decades is to deploy carbon capture and sequestration (CCS) in a mass scale [1,2]. The 
United States Department of Energy, and engineers, geologists, and scientists the world over have 
investigated the potential for reducing atmospheric carbon dioxide emissions through underground 
storage or geologic carbon sequestration. Although considerable research has been completed, no one 
has studied the potential for an integrated CCS system in Florida, USA. Florida is heavily dependent 
upon fossil fuels for its primary power generation. According to 2005 Environmental Protection 
Agency (EPA) data, there are 136 large and small power plants in Florida including those that use coal, 
oil, natural gas, and biomass as their source of fuel. Due to the large number of fossil fuel power plants 
in Florida, Florida also generates a significant mass of carbon dioxide emissions each year. In 2005, 
those emissions totaled almost 134,000,000 tonnes or 134 Mt [3]. In the future, new regulations are 
likely to require carbon dioxide emission reductions across the United States and in Florida. Florida 
utility companies are currently investigating feasible alternatives to reduce overall carbon dioxide 
emissions but most alternatives will be very expensive and difficult to implement over the short-term. 

While each alternative has advantages and disadvantages, forward thinking companies are likely to 
consider CCS due to its potential to sequester or store massive amounts of carbon dioxide and other 
greenhouse gases and its relatively low-cost attributes [4]. So what is CCS and how will it help reduce 
carbon liability? CCS is a technological innovation whereby carbon dioxide off-gas is captured, 
separated from other gases, concentrated, compressed, and then injected into underground repositories. 
Here the carbon dioxide is sequestered or stored for hundreds to thousands of years, effectively 
reducing the carbon footprint of the industrial emitter. In 2005 83% of Florida's electrical energy was 
produced by fossil fuels while in 2010 the percentage was almost 89% [5]. The continuing use of fossil 
fuels may depend upon finding suitable subsurface sequestration repositories in Florida and connecting 
them to a optimized network of pipelines and primary CO2 sources. 

According to the Intergovernmental Panel on Climate Change (IPCC), storage of CO2 in geologic 
formations includes four primary storage repository categories: saline aquifers, existing oil fields, 
depleted natural gas fields, and thin-nonmineable coal seams [6]. The capacity of each of these 
repository categories to sequester CO2 is an important planning variable to be considered during 
feasibility-level investigations of potential projects [1,7]. Deep saline aquifers appear to offer the 
highest potential capacity of the four primary options [8-11]. In addition, in Florida, saline aquifers are 
the most likely storage option. According to the United States Department of Energy [12], the 
estimated capacity of oil/gas fields is relatively small by comparison (e.g., 100 times less) and their 
geographic distribution is rather limited. A typical CCS saline aquifer storage project will undergo 
several operational changes over time with the injected CO2 ultimately becoming completely dissolved 
in the aquifer fluid. The various operation phases include site characterization, initial active injection, 
post-injection, and long-term monitoring. During the project lifecycle, there are significant changes in 
the state of injected CO2 with it starting as a free-phase, becoming residually-trapped, being dissolved, 
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and ultimately being precipitated as a mineral. The relative time scales for each process are different 
with residual trapping likely a decadal time scale, dissolution over hundreds of years, or more likely in 
saline waters, thousands of years and mineralization over even longer periods. During active 
operations, when liquid or supercritical CO2 is being injected into a repository, the CO2 will be highly 
mobile as a pure separate phase and concentrated aqueous phase [13]. CO2 is a highly compressible 
fluid compared to water and its density radically increases from 300 to 800 kg/m at pressure ranging 
from 10 to 25 MPa [14]. Since liquid or supercritical CO2 has a density less than the typical density of 
the saline repository fluid [15], it will be buoyant, tending to rise within the formation [16] until it 
intercepts a competent confining unit where it may spread laterally until it becomes trapped [17]. In 
some cases, depending upon formation dip, the supercritical CO2 may migrate updip along the 
confining unit. 

This paper describes ongoing research being conducted at the University of North Florida in 
Jacksonville, Florida, which is evaluating the feasibility of CCS projects within the state. The 
preliminary research results of this effort are presented in this paper and include an evaluation of an 
optimal Florida-wide pipeline transportation network connecting the 40 largest sources of CO2 
emissions with 5 alternate CO2 repositories or disposal areas. The 5 alternate repositories include a 
variety of saline aquifer storage zones located throughout Florida, including one offshore site located 
beneath the Gulf of Mexico, generally selected due to their favorable geologic properties. The 
optimized pipeline network is devised using "minimum cost transportation network model" 
methodology [18]. The use of this type of model for CO2 transmission in pipelines is relatively new 
and has not been investigated in Florida. Therefore, the first time application of the model 
methodology for this project in Florida is unique and provides a sound basis for more detailed 
evaluations in the future. 

2. Experimental Section 

2.1. Florida Primary Emission Sources 

The first task in developing an optimal CO2 pipeline transportation network for Florida, is to 
identify the location and magnitude of the largest sources of CO2 within the state. Florida has 
approximately 136 primary sources of CO2 inventoried by the EPA. For the initial model development 
effort, the 40 largest sources of CO2 were identified and summarized [3]. These 40 sources comprise 
over 90% of the 2005 total CO2 emissions for Florida. The 40 sources along with a map identification 
number, location in UTM 1983 (meters) horizontal grid coordinates, and the respective annual CO2 
emissions are listed in Table 1 . Each of the 40 sources is also shown on Figure 1 along with 5 potential 
CO2 repositories discussed later in this paper. 
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Table 1. 40 Largest sources of CO2 emissions in Florida [3]. 
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Florida Power Plant Location Map 




2.2. Florida Pipeline Transportation Model 

Following the identification of the major emission sources and magnitudes, the pipeline network 
model must be developed. The first step in this effort is to develop the pipeline cost model. Multiple 
investigators have developed and presented CO2 pipeline transportation models. Heddle et al. 
developed the initial MIT CO2 transport model [19]. Their model was a simple linear model that 
includes initial capital cost and annual operation and maintenance (O&M) costs. Equation (1) presents 
their model followed by an explanation of the variables. 

[Total Annual Cost = a xD xL x CF + O&M Cost] (1) 

where a is a constant estimated by MIT from literature values to be $33,853 (2003 dollars); 
D is the pipeline diameter in inches and is function of flow rate; 
L is the least-cost pipeline route length in miles; and, 
CF is a capital cost factor used to annualize the initial capital cost. 

McCoy developed a more sophisticated model that provides for regional cost differences as well as 
further resolution of cost factors such as pipe materials, labor, real estate, permitting, design, and 
construction management [20]. McCoy developed the cost model using regression analysis of natural 
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gas pipeline construction projects as published in the Oil and Gas Journal between 1994 and 2003 
McCoy model is presented in Equation (2) herein. 

[Total Annual Cost = (fi m x L a6m x D a7m x CF) + (fi L x L a6L x D a7L x CF) + 
(p RE x L a6RE x D a7RE x CF) + (fi MS x L a6MS x D a7MS x CF) + O&M Cost] 

where (3 m , (3l, (3re and (3 M s are cost coefficients for materials, labor, real estate, and miscellaneous 
(e.g., design, permitting, construction management) in 2004 dollars; L is the least-cost pipeline route 
length in kilometers; D is the pipeline diameter in inches and is a function of flow rate, inlet pressure, 
outlet pressure, and frictional losses; CF is a capital cost factor used to annualize the initial capital cost; 
a6 m , &6l, a6RE, and a6Ms are model pipeline length power exponents for materials, labor, real estate, 
and miscellaneous (discussed further below); and, a7 m , a7L, a7RE, and a7Ms are model pipeline 
diameter power exponents for materials, labor, real estate, and miscellaneous (discussed further below). 

In addition to MIT and McCoy, several others have developed recent cost models including 
Skovholt, Bakken & Von Streng Velken, and Zhang et al. [21-23]. Other recent papers have focused 
upon more sophisticated CCS infrastructure models [24-28]. These models simulate the full 
CCS value chain including carbon capture, compression, transport, and storage. For example, 
Mendelevitch et al. (2010) includes the "full decision path" for CCS infrastructure decisions in 
addition to consideration of a carbon tax or CO2 credits [24]. Kuby et al. (2011) discusses a CCS 
simulator that "determines the optimal quantity of CO2 to capture and optimize the various 
components of a CCS infrastructure network, given the price per tonne to emit CO2 into the 
atmosphere" [29]. All of the models are similar with capital costs dependent upon pipeline diameter 
and length. A key part of these models is the methodology to estimate the pipeline diameter. For most 
of these models, pipeline mass flow rate and diameter is a function of inlet pressure, outlet pressure, 
frictional losses, and topography (or elevation change). Since Florida is relatively flat, topographic 
differences are expected to be minor. In addition, if general transportation pressures are assumed to be 
consistent with industry practice, the development of a generalized method of estimating the required 
pipe diameter for a given mass flow rate is feasible. With the exception of Bakken & Von Streng 
Velken, the models that have been discussed in this paper so far have focused on landbased pipeline 
networks [22]. 

The IPCC presents pipeline transport costs over land and underwater versus CO2 mass flow rate [6]. 
This figure is reproduced herein for discussion purposes as Figure 2. As shown on the figure the cost 
relationship is highly nonlinear and the cost envelop is not constant. At a mass flow rate of 5 Mt per 
year, the cost envelope ranges from a lower bound of $2.10 for land construction and an upper bound 
of $4.50 for underwater construction. In addition, the costs for underwater pipeline are considerably 
more expensive per kilometer than land construction. In comparing the lower bound for land 
construction with the lower bound cost for underwater construction, underwater construction is 50 to 
75% more expensive. Bakken & Von Streng Velken present a cost model for a CCS project in Norway 
that is completely planned to be underwater [22]. In comparing these unit costs to the MIT model or 
McCoy model costs above, the unit cost is more than 2.5 times greater than equivalent land pipeline 
construction. In order to develop a new cost model for Florida, the percentages of land and underwater 
pipeline need to be calculated and a cost differential applied. This cost factor was assumed to be 
1.75 and is discussed further below. 



. The 

(2) 
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Figure 2. CO2 Pipeline transportation costs per tonne for 250 kilometer pipeline versus 
mass flow rate in Mt C0 2 per year (After IPCC (2005, Chapter 4, Figure 4.5)). 
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For this paper, as part of pipeline cost model development, the authors noted that all of the cost 
models require considerable calculations in order to estimate the required pipeline diameter. These 
calculations are certainly required for detailed pipeline design efforts but simpler estimates probably 
will suffice for planning or feasibility studies. The authors plotted pipeline diameter versus mass flow 
rate for several of the models discussed in order to develop a more generalized feasibility-level 
approach for estimating the pipeline diameter. Figure 3 presents the results of the various models. The 
authors fit published data to similar power models to ascertain a reasonable generalized model for the 
current research effort. The curve fitting statistics or R coefficients for the various power models were 
all greater than 0.99. For this research effort the authors chose to use a model between Skovholt and 
McCoy [20,21]. It is also interesting to note that the Zhang et at model returns practically the same 
values as the Heddle et a/./MIT "upper bound" diameter model [19,23]. The new proposed model 
provides a simplistic way to estimate the necessary pipeline diameter given solely an estimate of the 
required CO2 mass flow rate. The new model was developed for both English and SI units to permit 
practitioners from around the globe to use the model for preliminary planning purposes. Equations (3) 
and (4) provide the diameter estimation models in English and SI units, respectively. 
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Figure 3. CO2 Pipeline Diameter (inches) versus CO2 Mass Flow Rate Megatonne (Mt) 
CO2 per year for Multiple Published Cost Models (DIA = Diameter on Axis label). 
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[Pipeline Diameter, D = (11.25 x £ 0 - 3875 )] (3) 

where s is the CO2 mass flow rate in Megatonnes (Mt) per year; D is the pipeline diameter in inches. 

[Pipeline Diameter, D = (0.2857 x £ 0 - 3875 )] (4) 

where s is the CO2 mass flow rate in Mt per year; D is the pipeline diameter in meters. 

Other considerations for a new pipeline cost model are also important. First, the previous models 
were all developed in different years so that inflation adjustments are necessary to bring all models to 
2010 costs. Lewis provides construction cost factor data from 1996 to March 2010 for skilled labor, 
common labor, and materials [30]. From April 2004 to March 2010, the construction cost adjustment 
factors are 1.18 for materials, 1.25 for common labor, and 1.26 for skilled labor (e.g., for designers, 
permit specialists and construction managers). Real Estate costs have receded to close to 2005 costs 
such that a cost adjustment factor of only 1 .05 was assumed for the new model. In addition, pipelines 
may ultimately be required on land and underwater so a cost differential factor must be applied for 
underwater pipelines. As noted previously, the value of 1.75 was adopted for this study to convert a 
pipeline cost on land to one underwater. This value represents a reasonable average factor as 
determined from the literature. Lastly, the actual pipeline O&M cost needs to be determined. For 
detailed design, precise estimates of energy requirements, system maintenance, and inspection are 
required. For preliminary modeling purposes it is acceptable to use mean cost values from the 
literature as was done for this research effort. Using normalized unit O&M costs in $/tonne 
C02/kilometer from Zhang et al. and Bakken & Von Streng Velken, the authors developed a 
reasonable mean O&M cost of 0.0088 $/tonne C02/kilometer [22,23]. In reviewing the IPCC costs 
shown in Figure 2 above, one can calculate a similar mean O&M cost for comparison purposes [6]. 
Normalized O&M costs range from 0.005 $/tonne C0 2 /km for 30 Mt per year to 0.014 $/tonne 
C0 2 /km for 5 Mt per year. Mendelevitch et al. (2010) uses a value of 0.014 $/tonne C0 2 /km [24]. 
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Therefore, the adopted value of 0.0088 $/tonne CC^/km seems reasonable for a feasibility-level model 
especially since it is based upon more recent data than the original IPCC figure. With all of the 
required cost model elements assembled, the authors chose to adapt the McCoy model for use in 
Florida since it was thought to be the most complete and flexible [20]. Equation (5) provides the 
proposed new Florida C0 2 pipeline cost model. 

[Total Annual Cost = [[co m x p m x L a6m x (D x 39.38) a7m x CF] + [co L x p L x L a6L x (D x 39.38) a7L x CF] + 

[oo RE x p RE x L a6RE x (D x 39.38) a7RE x CF] + [m MS x p MS x L a6MS x (D x 39.38) a7MS x CF] + (5) 

[0.0088 x 62 x L] ] x a] 

where co m , co L , core and co M s are cost adjustment coefficients to convert April 2004 costs to March 2010 
costs and are co m = 1.18, col = 1.15, core = 1-05, and coms = 1-26; p m , Pl, Pre, and p M s are cost 
coefficients for materials, labor, real estate, and miscellaneous (e.g., design, permitting, construction 
management) in 2004 dollars and are p m = 1,534.62, p L = 30,690.22, p RE = 8,912.51, and 
Pms = 33,265.96; L is the least-cost pipeline route length in kilometers; D is the pipeline diameter in 
meters and is a function of flow rate (see Equation (4) above); CF is a capital cost factor of 0.067574 
assuming a 5% discount rate used to annualize the initial pipeline capital construction cost; s2 is CO2 
mass flow rate in tonnes per year; a is a factor to adjust costs for underwater construction, it is 1.75 for 
underwater projects and 1.0 for land pipeline projects; a6 m , a6L, a6RE, and a6MS are model pipeline 
length power exponents for materials, labor, real estate, and miscellaneous and are a6 m = 0.901, 
a6 L = 0.82, a6 RE = 1.049, and a6 M s = 0.783; and, a7 m , a7 L , a7 RE , and a7 MS are model pipeline diameter 
power exponents for materials, labor, real estate, and miscellaneous and are a7 m = 1.59, a7 L = 0.94, 
a7RE = 0.403, and a7jvis = 0.791. The new cost model for Florida is intended for use as a planning tool 
to be used in feasibility-level studies. It is applicable for use in Florida or other areas of similar 
flat topography. 

2.3. Geologic Repository or Disposal Zones 

With the pipeline cost model developed and the sources or supply nodes identified, the CCS 
repository or demand locations are identified next. The location was based upon the available geology, 
location of existing emission sources, and institutional concerns regarding possible CO2 releases [31]. 
Based upon the existing research, Florida has ample potential CCS repositories including depleted 
oil/gas fields, unminable coal seams, and deep, saline aquifers [32-38]. Of the four primary disposal 
alternatives, saline aquifers present the best opportunity to store large quantities of CO2 safely [4,12]. 
Building upon the existing research, this paper has chosen 5 separate saline aquifer CCS repository 
sites (see Figure 1) distributed throughout Florida and Southeast Georgia. Each of the 5 sites represents 
a portion of an identified CO2 disposal/repository site outlined in the "2010 Carbon Sequestration 
Atlas of the United States and Canada" [12]. The actual capacity of each of these prospective areas 
needs to be confirmed through further numerical modeling and subsurface testing [39,40]. Capacities 
assumed for the optimization modeling in this paper are on the low end of regional estimates. Each of 
these 5 sites is discussed herein. 

The Florida panhandle contains ample potential capacity for carbon sequestration within the Upper 
Cretaceous Zone, specifically the Tuscaloosa Formation. This formation is present in several Gulf 
Coast states and is estimated to have a "low" estimate capacity of at least 5 Gt according to the [12]. 
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For this study, the capacity within Florida (DAI on Figure 1) was estimated by the area- weighted share 
of the total estimated low capacity or 215 Mt. The top of the Tuscaloosa Formation is usually found 
below the shales of the overlying Eutaw Formation and in the area of Cedar Keys, Florida consists 
dominantly of red, light red, brown or mottled shales with interbedded sandstones [32,38]. Raymond 
and Copeland state in the Coastal Plain Province of Alabama, the Tuscaloosa Group comprises mainly 
fossilferous, nearshore, marine elastics [37]. In eastern-most Alabama, the formation is typically 
poorly sorted kaolinitic, arkosic sand and gravel interbedded yellowish-orange to reddish-green 
mottled kaolinitc clay. Thickness of the Tuscaloosa Formation ranges anywhere from 100 to 
400 meters [37]. USDOE describes the proposed storage reservoir at Southern Company Plant Daniel 
in Mississippi as a "massive sandstone that is a thick, regionally extensive, porous and permeable 
coastal to deltaic-marine sandstone at the base of the lower Tuscaloosa" [12]. According to the report, 
the lower Tuscaloosa in this area is overlain by a thick section of 90 to 140 meters of shales and 
mudrocks that were deposited as sea level rose during a marine transgression. 

DA2, shown on Figure 1, covers southeast Georgia. The proposed repository, the South 
Carolina-Georgia Basins, is located in saline aquifers that exist in the upper cretaceous layers of 
southeast Georgia, including the Atkinson Formation. The total preliminary "low" estimate capacity of 
this repository zone is 12.6 Gt [12]. For this study, the capacity (DA2 on Figure 1) was estimated by 
the area-weighted share of the total estimated low capacity or 4.98 Gt. The Atkinson formation is 
typically split into two separate units [34] . The separation of the units occurs roughly 990 meters from 
the surface [34]. The top unit tends to consist of shale which is micaceous, glauconitic and 
phosphoritic with the occasional non-uniform grain size sandstones. The lower unit of the Atkinson 
Formation, or the Basal Sand Unit, is a quartz sandstone, fine to coarse grained, poorly sorted, soft and 
typically glauconitic, phosphoritic and micaceous [34]. The Basal Sand Unit does contain sporadic 
stringers of shale, silt, and limestone. It is assumed that C0 2 disposal would primarily focus on the 
Basal Sand Unit. 

The repositories located in the peninsular region of Florida (DA4) and the offshore location (DA3) 
would inject CO2 into the Cedar Keys/Lawson Dolomite formations. USDOE estimates that the entire 
Cedar Keys/Lawson Dolomite formations capable of storing CO2 have a "low" estimate capacity of 
approximately 11 Gt [12]. For this study, the capacity of each area (DA3 and DA4 on Figure 1) was 
estimated by the area- weighted share of the total estimated low capacity or 1 Gt. According to Chen, 
the Cedar Keys Formation is widely spread across peninsular Florida and spreads into the 
panhandle [33]. In Brevard County, Florida, the top of the Cedar Keys Formation ranges from 
approximately 670 meters NGVD to 914 meters NGVD below land surface. The formation consists of 
dolomite and evaporates with a minor amount of limestone. Gypsum commonly fills pore spaces 
within the dolomite beds and occurs as thin irregular streaks or seams in the dolomite [33]. The 
Lawson Formation is generally found at the base of the Cedar Keys Formation. The Lawson is 
comprised mainly of pure, clean, very light brown and fine crystalline dolomite and/or chalky 
dolomitic limestone [33]. 

The disposal area indicated in south central Florida (DA5) has the potential to utilize multiple 
"stacked" secondary formations as a repository. USDOE estimates that the entire Cedar Keys/Lawson 
Dolomite formations capable of storing CO2 have a "low" capacity of approximately 11 Gt [12]. For 
this study, the capacity of the area (DA5 on Figure 1) was estimated by the area-weighted share of the 
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total estimated low capacity or 1 Gt. Including other stacked zones in the Paleocene and Eocene would 
contribute additional capacity for storing CO2, however, in order to provide a conservative estimate, 
this capacity was not included in the capacity estimate. In south Florida, zones of high transmissivity 
exist within the dolomite units of the Eocene, Paleocene, and Upper Cretaceous Formations [36]. The 
high transmissivity of the region is attributed to small to large cavities and caverns within the 
dolostone formations [35]. The cavities and caverns in the Upper Cretaceous Formation are typically 
found near the base of a massive dolostone facies in the upper part of the system and range in size 
from a baseball to a small car [36]. Data from nearby wells indicate the host rock is light to dark brown, 
anhedral dolestone. The same lithology holds true for the Paleocene and Eocene Cavity Zones [36]. 
The caverns tend to be surrounded by nearly impermeable material [35]. Utility companies local to the 
area have used the "Boulder Zone" for wastewater disposal in the past. Again, for the purpose of this 
study, these zones were not included in the storage estimates for area DA5. 

2.4. Optimization Model 

After the supply nodes and the demand nodes are identified and characterized, the transportation 
network is developed connecting supply nodes with demand nodes. Essentially, the least-cost 
transportation model can be developed using linear programming whereby a set of supply sources is 
linked to a set of demand locations with an optimal pipeline route where the unit costs of transportation 
between supply-demand pairs is known [18]. The primary model includes 40 possible CO2 supply 
sources each of which can be transported to any of 5 potential CCS demand locations or CCS 
repositories. Therefore, for the preliminary model there are 200 unknown CO2 "flows" for each year 
that must be solved using a set of linear equations. Each supply source can supply less than or equal to 
its annual CO2 emissions to the pipeline network. Each demand location can only accept less than or 
equal to its designated CCS capacity. The Florida pipeline transportation cost model is an adaptation of 
a logistics warehousing-type model developed by Cormier and Gunn and a CO2 pipeline optimization 
cost model developed by Bakken & Von Streng Velken [22,41]. This model is discussed below. The 
basic model equation and model constraints are included herein: 

[Minimize Y, Fijk x ^ siDjk = Total Cost] ^ 

where X is the annual CO2 pipeline transportation cost ($/tonne CO2) from CO2 supply node Si 
(from i = 1 to 40) to demand node or repository Dj (from j = 1 to 5) at Time Year k (from k = 1 to 
25 years) and Fyk is the CO2 flow through that pathway in tonnes CCVyear during Year k. 

[Subject to Constraint 1 Fijk < Capacity Dj] Summed from 1:200 each year (7) 
[Subject to Constraint 2 Fijk < Emission Supply from S t ] (8) 

[Subject to Constraint 3 Fijk > 0] (9) 

Using the new Florida cost model above, 200 annual CO2 pipeline unit transportation costs (X) 
were developed as part of the model development effort. These costs are shown on Table 2. 
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Table 2. Model CO2 pipeline unit transportation costs in Florida. 
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$7.79 
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1 1 
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C£ Q l 
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1Z 
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q> 1 / .4Z 


Co 1 A 

$y.i6 


$4.yu 


CO H£. 

$2. /6 


$Z.Z4 


13 


Port Everglades 


$19.49 


$14.31 


$8.33 


$4.23 


$1.37 


1 /i 
14 


Lansing Smith 
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$/.o4 
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15 


II I y' • 1 1 , 1 
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$14.00 
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C A C\£ 

$4.96 
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1 £ 

16 
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<£ a a T) 
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$7.25 


C c G£ 

$5.56 


c i /;o 

$1.62 


& A 1 (\ 
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1 / 
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<too 1£\ 
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cn o a 
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C/1 *7/1 

$4. /4 


C 1 A O 

$1.42 
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$4.56 


19 
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$10.52 


$6.02 
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20 
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$20.18 
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$2.55 


21 
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$12.05 
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$1.87 
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22 


P L Bartow 
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thn £L 1 
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23 
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$3.32 
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$6.57 
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24 
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$20.45 


$10.15 


$4.18 


$2.66 


25 
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Generating Station 
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26 
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$5.85 
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28 
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29 
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30 
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$7.86 


$2.40 
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Table 2. Cont. 
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36 


Suwannee River 


$13.92 


$8.98 


$25.06 


$13.67 


$20.38 


37 


Brandy Branch 


$18.71 


$6.90 


$24.85 


$11.92 


$19.36 


38 


Central Power & 


$30.05 


$14.89 


$14.36 


$6.86 


$13.76 


Lime 












39 


Putnam 


$21.48 


$10.22 


$20.70 


$9.18 


$17.20 


40 


Orlando Cogen LP 


$34.22 


$17.86 


$17.12 


$3.56 


$12.54 



Each of the unit costs was used to simulate the least-cost transportation network. In this first 
modeling effort, the least-cost transportation network was not constrained by geography, real estate 
limitations, institutional concerns (e.g. location of wetlands, parks, sensitive natural areas), or practical 
engineering considerations regarding pipeline right-of-way selection. Each possible emission 
source (SO was connected to each possible disposal repository (Dj) creating five possible pathways 
from each emission source. The optimization model calculated the least-cost alternative for each 
source-disposal pair. An example of this logic is shown on Figure 4 using the Crystal River plant as an 
illustrative emission source. Each pipeline route is broken into landward and underwater segments to 
demonstrate the importance of each part in the overall cost structure with underwater pipelines costing 
almost twice that of landward pipelines. The unit costs for each route range from $0.99 to 36.23 
per tonne CO2 with the arithmetic mean value at $10.49 per tonne per CO2. For the least-cost model as 
determined through the optimization effort, the total levelized cost for the network ranged from $2.29 
(during Year 1 for example) to $3.99 per tonne CO2 (during Year 25) which is in the range of $1 to $7 
per tonne CO2 referenced by McCoy [20]. As repositories are filled up beyond 25 years, the levelized 
cost could rise as high as $6.18 per tonne CO2. 
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3. Results and Discussion 

Following model development and testing, the preliminary optimal pipeline network was 
determined. Due to the high CO2 disposal area capacity at each proposed Florida repository location, 
the least-cost network model resulted in the use of all 5 repositories. Initially, in years 1 to 13, the 
proposed offshore repository was not utilized at all in the optimal network due to high unit costs for 
underwater pipeline installations. The underwater cost differential would have to be reduced by more 
than 30% for the offshore repository to be cost-effective during this timeframe. However, at the 
beginning of year 15, DA4 capacity is filled and CO2 flows initially directed to DA4 shift to DA2, 
DA3, and DA5. During years 1 to 13, the model indicated that approximately 7.62% of the modeled 
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emissions were transported to DAI, 14.05% to DA2, 58.63% to DA4, and 19.70% to DA5. The results 
for years 1 to 13 are shown in Table 3 and graphically presented on Figure 5. In year 14, as DA4 fills 
up, the model indicated that approximately 7.62% of the modeled emissions were transported to DAI, 
34.09% to DA2, 7.28% to DA3, 31.32% to DA4, and 19.70% to DA5. In years 15 to 22, the model 
indicated that approximately 7.62% of the modeled emissions were transported to DAI, 38.39% to 
DA2, 16.45% to DA3, and 37.54% to DA5. During year 23, DAI is filled to capacity. In year 23, as 
DAI fills up, the model indicated that approximately 3.04% of the modeled emissions were 
transported to DAI, 42.97% to DA2, 16.45% to DA3, and 37.54% to DA5. During the last 2 years of 
the model simulation, only three repositories are available and the remaining flows that went to DAI 
instead are directed to DA2. During the initial years of the simulation (years 1 to 13) none of the 5 
initial disposal areas is filled within the simulation period, repository capacity was not a model 
constraint so each of the 4 sites that were utilized received CO2 from emission sources closest to that 
site. For example, DA2 received all of its CO2 from large power plants in the Jacksonville, Florida area. 
Likewise, DA5 received all of its CO2 from large power plants in the Miami-Fort Lauderdale, Florida 
area. Additional research is warranted regarding maximum injection rates that can be sustained in each 
repository area as this may also affect the ultimate site feasibility. The optimal network included 5 out 
of the 5 simulated repository sites. If only one of the sites could be permitted and developed [42], it 
appears the most favorable site within Florida would be DA4, located south of Interstate 4 in Osceola 
County, Florida. This repository site will be further explored in future modeling to be completed as 
part of this research effort since under current capacity estimates; this site is filled within 14 years. One 
additional modeling scenario was generated that assumed that repositories 1, 4, and 5 were either filled 
or could not be permitted. Under these assumptions, almost 65% of the CO2 emissions would be 
transported to DA2 with the remainder going to DA3. Since a large flow is targeted to DA3, the 
offshore site, unit levelized costs would almost triple to $6.18 per tonne of C0 2 . In addition, future 
modeling efforts will evaluate an alternate pipeline network for Florida that will likely run along 
existing transportation right-of-ways (ROW) within Florida. This network will also combine smaller 
emission sites into "emission zones" where smaller contributing power plants will simply tap into 
larger diameter CO2 pressure mains emanating from the major emission loads. It is expected that this 
next generation optimization model will further reduce the overall projected cost for a Florida-wide 
pipeline network. Future modeling efforts may also need to include allowances for emission growth as 
Florida has experienced in the past. Florida's gross GHG emissions are rising faster than those of the 
nation as a whole (gross emissions exclude carbon sinks, such as forests). Florida's gross greenhouse 
gas emissions increased by about 35% from 1990 to 2005, while national emissions rose by 16% from 
1990 to 2005 [43]. The growth in Florida's emissions from 1990 to 2005 is primarily associated with 
the electricity consumption and transportation sectors [43]. The effect of alternative energy and nuclear 
projects under development in Florida may also impact the analysis by actually reducing the overall 
emissions. Future modeling efforts should also include a factor that weights the emission sources by 
their probability of using CCS. Also, future models should also evaluate the actual surface "footprint" 
of the CO2 in order to refine the estimated repository capacity [44]. 
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Table 3. Preliminary optimal CO2 pipeline network for Florida. 



Map ID 


Plant/Facility Name 


Annual C0 2 
Emission (Mt) 


Optimum Disposal Area 
in Florida (DA 1 to 5) 
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9.13 
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3 
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8 


Sanford 
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5.12 
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10 
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4.96 
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11 
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3.91 
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12 
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3.79 


DA5 
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4. Conclusions 

In conclusion, this paper discusses the development of and preliminary results from, a model of a 
Florida-wide CO2 pipeline transportation network. The top 40 primary CO2 emission sources in 
Florida were linked to 5 hypothetical repositories spread throughout the state including one offshore 
site within the Gulf of Mexico. The optimization model determined that an optimal network of 5 sites 
results in the lowest preliminary network cost which is estimated to cost approximately $289,000,000 
to $503,000,000 per year for 25 years. During years 1 to 13 when all repositories are available for 
disposal, the total annual levelized cost is $2.26 per tonne CO2. As repository capacity fills up, this 
levelized cost rises to $2.58 per tonne in year 14; $3.52 per tonne in years 15 to 22; $3.76 per tonne in 
year 23; and, $3.99 per tonne in years 24 to 25. Beyond the 25-year simulation period, it is estimated 
that only DA2 and DA3 would have remaining capacity to store CO2 and the levelized cost would 
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increase further to $6.18 per tonne CO2. This preliminary modeling effort suggests that a fully- 
connected pipeline network would not be the lowest cost CCS plan as long as all five repository 
locations are available. 

Acknowledgements 

This research was supported by the University of North Florida (UNF) College of Computing, 
Engineering, and Construction Management as well as the Academic Affairs Unit from UNF. The 
researchers would also like to thank Greg Powell from Golder Associates, Inc. for his detailed review 
of the draft manuscript. 

References and Notes 

1. Bradshaw, J.; Bachu, S.; Bonijoly, D.; Burruss, R.; Holloway, S.; Christensen, N.P.; 
Mathiassen, O.M. C0 2 storage capacity estimation: Issues and development of standards. Int. J. 
Greenh. Gas Control 2007, 1, 62-68. 

2. The Future of Coal: Options for a Carbon Constrained World; Massachusetts Institute of 
Technology: Cambridge, MA, USA, 2007. Available online: http://mit.edu/coal/ (accessed on 
15 September 2010). 

3. Environmental Protection Agency. Year 2005 eGRID Plant, Boiler, and Generator Data Files; 
EPA Clean Energy Resources: Washington, DC, USA, 2007. Available online: http://epa.gov/ 
cleanenergy/energy-resources/egrid/index.html (accessed on 18 August 2010). 

4. United States Department of Energy (DOE). Methodology for Development of Geologic Storage 
Estimates for Carbon Dioxide, 2008; National Energy Technology Laboratory: Pittsburgh, PA, 
USA, 2008; pp. 115-132. 

5. Energy Information Administration (EIA). Florida State Energy Profile Data 2009 and 2010. 
Available online: http://www. eia.gov/state/state_energy_profiles. cfm?sid=FL (accessed on 1 
December 2010). 

6. Benson, S.; Cook, P. Underground Geological Storage, Special Report on Carbon Dioxide 
Capture and Storage; The Intergovernmental Panel on Climate Change (IPCC), United Nations: 
New York, NY, USA, 2005; pp. 1-265. 

7. Koide, H.G.; Tazaki, Y.; Noguchi, Y.; Nakayama, S.; Iijima, M.; Ito, K. Subterranean 
containment and long-term storage of carbon dioxide in unused aquifers and in depleted natural 
gas reservoirs. Energy Conv. Manage. 1992, 33, 619-626. 

8. Bachu, S.; Gunter, W.D.; Perkins, E.H. Aquifer disposal of CO2: Hydrodynamic and mineral 
trapping. Energy Conv. Manage. 1994, 35, 269-279. 

9. van der Meer, L.G.H. The CO2 storage efficiency of aquifers. Energy Conv. Manage. 1995, 36, 
513-518. 



Int. J. Environ. Res. Public Health 2011, 8 



973 



10. Obdam, A.; van der Meer, L.G.H.; May, F.; Kervevan, C; Bech, N.; Wildenborg, A. Effective 
CO2 Storage Capacity in Aquifers, Gas Fields, Oil Fields and Coal Fields. In Proceedings of the 
6th International Conference on Greenhouse Gas Control Technologies (GHGT-6), 1-4 October 
2002, Kyoto, Japan; Gale, J., Kaya, Y., Eds.; Pergamon: Oxford, UK, 2003; pp. 339-344. 

11. Herzog, H. Carbon Dioxide Capture and Stoage. In Economics and Politics of Climate Change; 
Helm D., Hepburn C, Eds.; Oxford University Press: Oxford, UK, 2009; pp. 263-283. 

12. United States Department of Energy (DOE). Carbon Sequestration Atlas III of the United States 
and Canada; Southeast Regional Carbon Sequestration Partnership, National Energy Technology 
Laboratory: Pittsburg, PA, USA, 2010; pp. 84-95. 

13. Bachu, S.; Adams, J.J. Sequestration of CO2 in geological media in response to climate change: 
Capacity of deep saline aquifers to sequester CO2 in solution. Energy Conv. Manage. 2003, 44, 
3151-3175. 

14. Han, W.S.; McPherson, B.J. Optimizing geologic CO2 sequestration by injection in deep saline 
formations below oil reservoirs. Energy Conv. Manage. 2009, 50, 2570-2583. 

15. Sharqawy, M.H.; Lienhard, J.H.; Zubair, S.M. Thermophysical properties of seawater: A review 
of existing correlations and data. Desal. Water Treat. 2010, 16, 354-380. 

16. Carbon Dioxide Thermophysical Property Calculator; Massachusetts Institute of Technology: 
Cambridge, MA, USA, 2010. Available online: http://sequestration.mit.edu/tools/index.html 
(accessed on 15 September 2010). 

17. Flett, M.A.; Gurton, R.M.; Taggart, I.J. Heterogeneous Saline Formations: Long-Term Benefits 
for Geo-Sequestration of Greenhouse Gases. In Proceedings of the 7th International Conference 
on Greenhouse Gas Control Technologies, Vancouver, Canada, 5-9 September 2004; Rubin, E.S., 
Keith, D.W., Gilboy, C.F., Morris, T., Thambimuthu, K, Eds.; Elsevier: New York, NY, USA, 
2005; pp. 501-509. 

18. Powell, S.; Baker, K. Management Science. In The Art of Modeling with Spreadsheets, 3rd ed.; 
John Wiley & Sons: Hoboken, NJ, USA, 2009; pp. 281-309. 

19. Heddle, G.; Herzog, H.; Klett, M. The Economics of C0 2 Storage; Report # 2003-003; MIT Press: 
Cambridge, MA, USA, 2003. 

20. McCoy, S. The Economics of CO 2 Transport by Pipeline and Storage in Saline Aquifers and Oil 
Reservoirs; Ph.D. Thesis; Carnegie Mellon University Press: Pittsburg, PA, USA, 2008; 
pp. 1-247. 

21. Skovholt, O. CO2 transportation system. Energy Conv. Manage. 1993, 34, 1095-1103. 

22. Bakken, B.H.; von Streng Velken, I. Linear Models for Optimization of Infrastructure for CO2 
Capture and Storage. IEEE Trans. Energy Convers. 2008, 23, 824-833. 

23. Zhang, Z.X.; Wang, G.X.; Massarotto, P.; Rudolph, V. Optimization of pipeline transport for CO2 
sequestration. Energy Conv. Manage. 2006, 47, 702-715. 

24. Mendelevitch, R.; Herold, J.; Oei, P.-Y.; Tissen, A. CO2 Highways for Europe: Modelling a 
Carbon Capture, Transport and Storage Infrastructure for Europe; Centre for European Policy 
Studies: Brussels, Belgium, 2010; p. 23. 



Int. J. Environ. Res. Public Health 2011, 8 



974 



25 van den Broek, M.; Brederode, E.; Ramirez, A.; Kramers, L.; van der Kuip, M.; Wildenborg, T.; 
Turkenburg, W.; Faaij, A. Designing a cost-effective CO2 storage infrastructure using a GIS based 
linear optimization energy model. Environ. Modell. Softw. 2010, 25, 1754-1768. 

26. Middleton, R.S.; Bielicki, J.M. A scalable infrastructure model for carbon capture and storage: 
SimCCS. Energy Policy 2009, 37, 1052-1060. 

27. Middleton, R.S.; Bielicki, J.M. A comprehensive carbon capture and storage infrastructure model. 
Energy Procedia 2009, 1, 1691-1698. 

28. Keating, G.N.; Middleton, R.S.; Stauffer, P.H.; Viswanathan, H.S.; Letellier, B.C.; Pasqualini, D.; 
Pawar, R.J.; Wolfsberg, A.V. Mesoscale. Carbon sequestration site screening and CCS 
infrastructure analysis. Environ. Sci. Technol. 2011, 45, 215-222. 

29. Kuby, M.; Bielicki, J.M.; Middleton, R.S. Optimal spatial deployment of carbon dioxide 
capture and storage given a price on carbon dioxide. Int. Reg. Sci. Rev. 2011, 
doi: 10.1177/0160017610397191. 

30. Lewis, S. 1Q Cost Report. Eng. News Rec. 2010, 264, 58. Available online: 
http://enr.construction.com/economics/historical_indices/ (accessed on 10 October 2010). 

31. Lewicki, J.L.; Birkholzer, J.; Tsang, C. Natural and industrial analogues for leakage of CO2 from 
storage reservoirs: Identification of features, events, and processes and lessons learned. Environ. 
Geol. 2007, 52, 457-467. 

32. Cole, W.S. Stratigraphic and Paleontologic Studies of Wells in Florida — No. 2. Florida 
Geological Survey; Geological Bulletin No. 20; Florida Geological Survey: Tallahassee, FL, USA, 
1942. Available online: http://ufdc.ufl.edu/fgs (accessed on 15 November 2010). 

33. Chen, C.S. The Regional Lithostratigraphic Analysis of Paleocene and Eocene Rocks of Florida. 
Florida Geological Survey; Geological Bulletin No. 45; Florida Geological Survey: Tallahassee, 
FL, USA, 1965. Available online: http://ufdc.ufl.edu/fgs (accessed on 15 November 2010). 

34. Babcock, C. Geology of the Upper Cretaceous Clastic Section Northern Peninsular Florida. 
Florida Geological Survey; Information Circular No. 60; Florida Geological Survey: Tallahassee, 
FL, USA, 1969. Available online: http://ufdc.ufl.edu/fgs (accessed on 15 November 2010). 

35. Vernon, R.O. The Beneficial Uses of Zones of High Transmissivity in the Florida Subsurface for 
Water Storage and Waste Disposal; Florida Geological Survey Information Circular No. 70; 
Florida Geological Survey: Tallahassee, FL, USA, 1970. Available online: http://ufdc.ufl.edu/fgs 
(accessed on 15 November 2010). 

36. Puri, H.S.; Winston, G.O. Geologic Framework of the High Transmissivity Zones in South 
Florida; Florida Geological Survey Special Publication No. 20; Florida Geological Survey: 
Tallahassee, FL, USA, 1974. Available online: http://ufdc.ufl.edu/UF00000159/00001/ 
56x?td=puri (accessed on 15 November 2010). 

37. Raymond, D.E.; Copeland, C.W. Alabama Stratigraphy; Geological Survey of Alabama Circular 
No. 140; Florida Geological Survey: Tallahassee, FL, USA, 1988. Available online: 
http://www.gsa.state.al.us/online_pubs.aspx (accessed on 15 November 2010). 

38. Rupert, F.R. Geology of Gulf County, Florida; Florida Geological Survey Bulletin No. 63; Florida 
Geological Survey: Tallahassee, FL, USA, 1991. Available online: http://ufdc.ufl.edu/fgs 
(accessed on 15 November 2010). 



Int. J. Environ. Res. Public Health 2011, 8 



975 



39. Yamamoto, H.; Zhang, K.; Karasaki, K.; Marui, A.; Uehara, H.; Nishikawa, N. Numerical 
investigation concerning the impact of CO2 geologic storage on regional groundwater flow. 
Int. J. Greenh. Gas Control 2009, 3, 586-599. 

40. Okwen, R.T.; Stewart, M.T.; Cunningham, J. Analytical solution for estimating storage efficiency 
of geologic sequestration of C0 2 . Int. J. Greenh. Gas Control 2010, 4, 102-107. 

41. Cormier, G.; Gunn, EA. Modelling and analysis for capacity expansion planning in warehousing. 
J. Oper. Res. Soc. 1999, 50, 52-60. 

42. Bayer, M.R.; Kobelski, B. Underground injection control program: Proposed regulations for 
underground injection of carbon dioxide for geologic sequestration. Ground Water Monit. 
Remedial 2008, 28, 42-44. 

43. Center for Climate Strategies. Final Florida Greenhouse Gas Inventory and Reference Case 
Projections 1990-2025; Technical Report for the Florida Department of Environmental Protection; 
Center for Climate Strategies: Washington, DC, USA, 2008; p. 120. 

44. Brennan, S.T.; Burruss, R.C. Specific storage volumes: A useful tool for CO2 storage capacity 
assessment. Nat. Resour. Res. 2006, 15, 165-182. 



© 2011 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.0rg/licenses/by/3.O/). 



